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ABSTRACT 


The application of a thyristor (a four-layer P1-N1-P2-N2 semiconductor struc- 
ture) as an optical detector is explored. Based on laboratory experiments which demon- 
strated that this device produces a pulse-mode output to incident light, the thyristor is in- 
vestigated by comparing the existing theory of static forward-biased operation to simula- 
tion results obtained using ATLAS by Silvaco, Inc. The results include identification of 
the holding point on the IV curve by simulating the junction potential across each junc- 
tion as a function of current, and demonstration that impact ionization is not a critical fac- 
tor in thyristor operation. A series of simulations were performed which show that the 
thyristor can be optimized for use as a detector by decreasing the emitter efficiencies by 
decreasing the doping in the P1 and N2 layers, or by increasing the doping in the P2 
layer; the switching voltage can be controlled by selecting the doping and thickness of 
the N1 layer. A detector device was designed to allow further testing of the thyristor de- 
tector using the ABN CMOS process from AMI Semiconductor via the MOSIS service. 


The design of this device is discussed and simulated IV curves are presented. 
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EXECUTIVE SUMMARY 


Existing semiconductor photodetectors produce a direct current output in response 
to incident light. The magnitude of the current output from the device is related to the in- 
tensity of the incident light. Various noise mechanisms exist in photodetectors which 
cause ransom variations in the output current magnitude. The work contained in this the- 
sis explores a novel photodetector mechanism which uses a thyristor (a four-layer P1-N1- 
P2-N2 semiconductor structure) in a simple RC circuit as an alternative to the photodi- 
ode. The demonstrated response of this circuit is a pulse-mode signal which varies in fre- 
quency with the intensity of the incident light. 

Thyristors have been used for many decades, with the typical application being a 
switch. Specially designed thyristors, called LASCRs, have been designed to be turned 
on using photo-generated carriers. The detector application of the thyristor depends, how- 
ever, on preventing the thyristor from completely turning-on. Much of the existing litera- 
ture on thyristor theory is intended to better understand how to enhance the turn-on proc- 
ess. A focus of this work is to identify those features which can be used to prevent com- 
plete turn-on of the thyristor, thereby enhancing its suitability as a detector. 

A simplified thyristor structure was used for simulations performed by ATLAS 
from Silvaco, Inc, to investigate the behavior of the device at different operating points 
under forward biasing. Two key operating points, the switching point where the thyristor 
begins turn-on and the holding point where the turn-on process is completed, are defined 
using simulation results of the electric field and potential at the pn junctions in the device. 
The switching point is defined as the point where the built-in electric field at the middle 
junction reaches its maximum strength. This, of course, corresponds with the maximum 
applied potential at this junction. The holding point is defined as the point where the ap- 
plied potential is reduced to zero. In other words, the junction potential is equal to the 
built-in potential of the junction. 

Impact ionization is often described as the current multiplication mechanism 
which gives the thyristor its characteristic shape. In this thesis, simulations and analysis 


of the electric field in the device show that impact ionization is not a critical factor. 
xili 


Simulations show that impact ionization does not occur in some thyristors and does not 
substantially change the IV curves of devices where it does occur. 

A series of simulations were performed to explore the effect of variations to the 
doping concentrations and thicknesses of the thyristor layers. The effort was to identify 
those relationships which enhance the operation of the thyristor as a detector. It is critical 
for a thyristor detector not to complete the turn-on process, which is to say that the cur- 
rent must be kept below the holding current. Simulations showed that the holding current 
for a given thyristor structure can be increased by lowering the emitter efficiency of the 
two outer pn junctions of the thyristor. This can be done by decreasing the doping con- 
centration of the two emitter layers, P1 and N2, or by increasing the doping concentration 
of the P2 base layer. 

Finally, a practical thyristor detector device was designed for fabrication using a 
CMOS process with a bipolar option. The design and layout for the device are discussed 
and simulations of the structure are presented. This device, while not an optimal detector 
structure, will permit controlled experiments to further develop the thyristor detector con- 


cept. 


XIV 


I. INTRODUCTION 


The purpose of this thesis was to investigate the static operation of thyristors, spe- 
cifically with the intent to use them as optical detectors. The goal that this thesis supports 
is a novel silicon photodetector structure based on the thyristor. Although thyristors and 
other semiconductor photodetectors have existed for many decades, the application pur- 
sued herein represents a departure from the standard direct current output of existing 
semiconductor photodetector devices. The thyristor, when placed in an appropriate cir- 
cuit, generates a pulsed output that varies in frequency with changes to the incident light 


intensity. 


This chapter gives a brief introduction to the physics of photon absorption, par- 
ticularly as it is occurs in photovoltaic detectors. This is followed by a synopsis of the 
experimental results which motivated this work. The second chapter of the thesis ex- 
plores the structure and physical operation of thyristors. The basic device is described 
along with an introduction to the existing static theory of operation for the device. The 
chapter concludes with some results obtained by computer simulations which are relevant 
to thyristors in general. The third chapter discusses the design and computer simulation 
of a detector device created specifically for use as an optical detector which will be fabri- 
cated using a commercially-available manufacturing process. The final chapter summa- 


rizes the results of the thesis and discuss possible areas of future research. 


A. PHOTON DETECTION USING SEMICONDUCTORS 

A photodetector can be defined as a class of optical detector which detects pho- 
tons by a measurable interaction between incident photons with the material, resulting in 
a change in the physical state of the material [1]. The physical property which changes in 
semiconductors is the free carrier concentration, provided that the photons have sufficient 
kinetic energy. The energy of a photon Epjoio, iS inversely related to its wavelength 2 by 


he 
Ee i = A (1.1) 


where h is Planck’s constant and c is the speed of light. When Epjoion exceeds the energy 
band gap E, of the semiconductor material, the photon can be absorbed. Absorption cre- 


ates an electron-hole pair of free carriers which can be used to conduct electric current. 


Semiconductor material can be used to make two varieties of photodetectors. The 
first is the photoconductor which responds to light of sufficient energy by reducing the 
resistance to current flow through additional current carriers. The change in resistance is 
detected by an external biasing circuit. The second, and more common, type is the photo- 
voltaic detector. This uses an internal electric field to sweep the photon-created carriers 
away before they can recombine. This is the basis for devices such as photodiodes and 
solar cells which uses the internal electric field found at unbiased or reverse-biased pn 
junctions. An important limitation to photovoltaic detectors is that only those carriers 
created in a region where an electric field is present are usable. Carriers created away 


from the electric field will quickly recombine and be of no further use. [1] 


As will be seen in the following chapter, a thyristor is nothing more than a pnpn 
structure which consists of three pn junctions in series. Under forward bias of the thyris- 
tor, one of those junctions is actually reverse biased. This extends the width of the elec- 
tric field, making that junction more useable as a photovoltaic detector. In this respect, 
the thyristor is then nothing more than a more complicated structure that performs the 
same function as a photodiode. However, the nature of the response of the thyristor is en- 


tirely different than a photodiode, as will be discussed in the following section. 


B. SUMMARY OF PREVIOUS EXPERIMENTAL RESULTS 

This section summarizes laboratory experiments by Karunasiri [2] which have 
demonstrated that a specially prepared thyristor placed in a simple RC circuit produces a 
pulsed output when illuminated with light. The thyristor used in this experiment was a 
Motorola MBS4993 which was a relatively small, commercially available device. It was 
not designed for this application, but rather as a bidirectional switch. To make the device 
useable for this application, the protective packaging was etched away exposing the ac- 
tual semiconductor device to incident light. This is not an ideal device on which to con- 


duct further tests. One of the objectives of this thesis is to design a device which will 
2 


provide a better testing platform. A general schematic diagram of the experiment circuit 


is given in Figure | [2]. 


Light 
Pulse stream 
inal 






Thyristor 
C R Output 


Figure 1. Circuit Diagram (After [2].). 


The pulsed output of this circuit is what makes the thyristor detector an interesting 
concept. The response of a photodiode is a direct current signal which changes amplitude 
with the intensity of the incident light. The pulse mode output, however, changes fre- 


quency with the intensity of incident light. 


An example of the pulse mode output is given in Figure 2. Notice that the fre- 


quency of pulses increases as the light power on the thyristor is increased. 
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Figure 2. Pulse Mode Output Variation with Incident Light Intensity (From [2].). 
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There are two aspects of the pulse mode output which make this device worthy of 
further investigation. The first is related to the potential noise advantage that the pulse 
mode signal may have over a direct signal. The second is the potential application of this 


detector. These will be discussed briefly below. 


Noise in photodiodes manifests itself in variations in the output current while the 
incident light is constant. Presumably the same noise sources will be at work in the thy- 
ristor photodetector, but the variations will result in changes to the pulse interval. Two of 
the more significant noise mechanisms in voltaic photodetectors are shot noise and gen- 
eration-recombination noise. These both involve the statistical variations in the time be- 
tween events, either between electrons passing across a potential barrier in the case of 
shot noise or between generation-recombination events in the case of generation- 
recombination noise, which can be described by Poisson statistics involving extremely 


large numbers of events. [1] 


The Poisson distribution approaches a Gaussian distribution around some mean 
when the number of events is large, as is usually the case with photons and electrons. 
The Gaussian distribution is then symmetrical about the mean. This means that the inter- 
val between individual events may vary, but the total interval measured over a large 
number of events will be approximately equal to the number of events times the mean in- 
terval. Thus, if the pulse interval is large compared to the interval between electron or 
photon events, then the statistical variation in event intervals will not significantly impact 


the pulse interval. [1] 


Since noise is not likely to cause significant pulse interval variation shown in Fig- 
ure 2, another explanation needs to be sought. That work can begin after a better testing 
platform is fabricated. At this point, it is not thought that the variation represents a prob- 
lem for use of the thyristor as a detector. This again assumes that the interval distribution 
is symmetrical about some mean. Then the variation can be mitigated by averaging the 


interval over a series of pulses. 


The second interesting characteristic is the pulse train’s resemblance to the wave- 


form used by the human body to transmit electrical signals over nerves. Specifically, the 
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output of the eye to light stimulation can be described as a “digital spike train” transmit- 
ted to the brain [3]. This points to one of the possible future applications of a thyristor- 


based photodetector, an array of photodetectors to act as a retinal prosthesis. 


This chapter has introduced photodetection using semiconductor materials in gen- 
eral, with a specific emphasis on the use of a thyristor to obtain a pulse-mode detector 
output. The following chapter introduces the structure of the thyristor and static thyristor 
theory. Simulations were performed to better describe the physical state of the thyristor 
during forward-biased operation. Finally, the design for a device to be fabricated to per- 


mit further testing is introduced. 
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Il. THYRISTORS 


This chapter presents an introduction to the structure and operation of the four- 
layer semiconductor device known as the thyristor. A brief overview of thyristor struc- 
ture and static operation is followed by an introduction to static thyristor theory. The 
chapter concludes with comprehensive static simulations of a thyristor in an effort to bet- 


ter understand thyristor operation. 


A. OVERVIEW 

1. Structure 

The name thyristor is a general term which applies to a large family of four-layer 
p-n-p-n semiconductor devices. A diagram of the basic structure of a thyristor is shown 


in Figure 3. 


Gate (s) 





Figure 3. _Thyristor Structure. 


Currents and voltages described in this thesis will follow the conventions shown 
in Figure 3. Forward biasing implies that V4x is positive, while positive anode current ly 
is defined as positive when flowing into the P1 layer via the anode contact. The literature 


uses a wide range of names for the different layers and parts of a thyristor. This thesis 


will adhere to the names used in Figure 4 and described below. The anode is the electri- 
cal contact to the P1 layer and the cathode is the electrical contact to the N2 layer. Gate 
contacts may be used to either, both, or neither, of the N1 and P2 layers. The junctions 
formed by the layers are labeled as follows: J1 is the junction between P1 and N1, J2 is 


the junction between N1 and P2, and J3 is the junction between P2 and N2. 


The number of gate contacts is one way of sorting thyristors. For example, a thy- 
ristor without gate contacts, that is with only anode and cathode contacts, is known as a 
Shockley diode; a thyristor with one gate contact to the P2 layer is known as a silicon- 
controlled rectifier (SCR); and a thyristor with contacts to both the N1 and P2 layers is 
known as a silicon-controlled switch (SCS). The devices described in this chapter will be 


Shockley diodes unless clearly specified otherwise. 


The relative doping levels of the layers and the widths of the two internal layers 
are important to thyristor operation. Some general guidelines for these parameters are 
given by Pierret [4]; the Pl and N2 layers should be highly doped relative to the N1 and 
P2 layers, and the thicknesses of the N1 and P2 layers should be on the order of a minor- 
ity carrier diffusion length in each layer. A representative doping profile of a thyristor is 
shown in Figure 4 where the vertical axis is a logarithmic scale. This plot format is used 
frequently in this thesis to show various concentration quantities. In each case the verti- 


cal scale is the logarithm of the actual concentration value. 
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Figure 4. Representative Thyristor Doping Profile (After [4].). 
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2. I-V Characteristic 
Thyristors have a characteristic current vs. voltage, or IV, curve, which is repre- 
sented in Figure 5. Again, the literature is inconsistent in the labeling of regions and 


points on the curve. This thesis will consistently use the labels shown in Figure 5. 


Vax 





Figure 5. Typical Current-Voltage Curve of a Thyristor (After [5].). 


The IV curve can be broken into several operating modes. These operating modes 


are described by considering the bias voltage V4x and the anode current /4 within the for- 


ward and reverse bias regions. Beginning with no bias, as reverse bias (V,, <0) is ap- 


plied, the thyristor initially allows only a small reverse current to flow. This is the re- 
verse blocking mode. Once the bias becomes more negative than the reverse breakdown 
voltage Ve the thyristor enters the reverse breakdown region. This is characterized by 
rapidly increasing negative anode current with increasingly negative bias. The reverse 


bias operation of a thyristor behaves similar to a standard pn junction under reverse bias. 


The forward bias region is more complicated. Because of the multiple current 
values, for most voltage points it is easiest to understand this region by considering the 
anode current rather than the anode-to-cathode voltage. Beginning from zero, increasing 


the bias voltage initially produces only a very slowly increasing anode current. When the 
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anode current is less than the switching current Js, the operating mode is called the for- 
ward blocking mode. It is characterized by very high resistance to current. The switch- 
ing voltage Vs is the maximum voltage across the device. When the current exceeds Is 
the device enters the negative resistance mode. As the name indicates, this mode is char- 
acterized by negative resistance where increasing anode current results in anode-to- 
cathode voltage lower than Vs. Once the anode current exceeds the holding current Jz, 
the device is considered to be on. This is the forward conducting mode and it is charac- 
terized by extremely low resistance. The holding voltage V7; is the local minimum con- 


current with J;,. [5] 


It is important to emphasize the distinction between the switching point and the 
holding point. Switching is not synonymous with turning on. The device begins turning 
on at the switching point, but is not fully on until the holding point is reached. Thyristors 
are usually designed to enhance the turning-on process and do not normally operate be- 
tween the switching and holding points. Returning the device to the forward blocking 
mode is much more involved from the forward conducting mode than from the negative 
resistance mode. 

3. Switching Mechanisms 

There are several ways to induce switching in a thyristor. The most direct is in- 
creasing the bias voltage until the anode current exceeds the switching current. This is 
how a Shockley diode is switched, but it also applies to other thyristors. As long as the 
biasing source is not current-limited below the holding current, the device will switch on. 


This mechanism will be used in many simulations throughout this chapter. 


Thyristors with a gate contact can be switched using an appropriate gate current 
which is much smaller than the anode current. This method is typically used in the clas- 
sic application of the thyristor as a switch. As an example, consider an SCR connected as 


shown in Figure 6. Recall that an SCR is a thyristor with a gate contact to the P2 layer. 
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Figure 6. Example SCR Circuit. 


Increasing the gate current results in a decrease in the switching voltage Vs 
as shown in Figure 7 where J,, < /,, <1,,. Since the gate current is orders of 
magnitude smaller than the anode current [4], the thyristor can be used to control 
very large anode currents and voltages with relatively small gate current. SCR’s 
can be designed to block many thousands of volts in the forward blocking region 


and to conduct many thousands of amperes in the forward conducting mode [5]. 
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Figure 7. | Reduction in Vs with increasing J. 
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Thyristors can also be switched by using incident light. Increasing incident inten- 
sity acts much like increasing gate current to lower the switching point, except the carri- 
ers are much more quickly swept through the depletion region because they do not have 
to first travel from the gate. Furthermore, photon absorption generates two carriers, both 


of which contribute to anode current. [6] 


Light activation of thyristor switches allows complete electrical isolation of the 
switch from the controlling circuit. Devices designed to behave as such are called Light- 
Activated SCRs (LASCR). LASCRs have a switching point that varies with incident 
light intensity, as described above. It must be emphasized that LASCRs are not suitable 
for use as optical detectors since they are designed to turn on when sufficient light energy 
is incident on the device. As mentioned above, turning off a thyristor operating above the 
holding current level is not a simple process. At best, a LASCR employed as a sensor 
will be able to indicate that at one point in time light was incident on the device. It will 
not be able to indicate when that may have occurred or even if the light is still incident. 
As such, it is not of much use as a detector. This points to a design objective for a thyris- 


tor used as photodetector; avoid features which enhance the turn-on process. 


B. THYRISTOR OPERATIONAL THEORY 

This section introduces thyristor operational theory. The focus of this section is 
limited to the forward-biased region since that is the region of interest in optical detection 
with thyristors. Coverage of reverse-biased theory can be found in Blicher [7] or Ghan- 
dhi [8]. Further, this thesis is interested only in the operation of the thyristor below the 
holding point. Issues related to completing the turn-on process and turning the device off 
are also not covered. Again, Blicker and Ghandhi provide good coverage of these topics. 

1. Introduction to Thyristor Theory 

The development of the thyristor did not lag far behind the development of the bi- 
polar junction transistor. Shockley described what became known as the thyristor in his 
seminal work on bipolar devices in 1950. He described the four-layer p-n-p-n structure 
as a high current-gain variation of the p-n-p transistor, labeling the additional N2 layer as 


the current multiplication layer [9]. 
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In trying to understand the physical reason for the high gain observed in the high- 
gain transistor, which occurs in the negative resistance operating region for a thyristor, 
Shockley effectively ruled out secondary generation, or avalanche multiplication, of car- 
riers. He asserted that the high gain occurs even when electric field strength in the deple- 
tion regions is insufficient to induce secondary generation. Shockley concluded his in- 
vestigation by stating that the likely cause of the high gain is what he calls the p-n hook 
theory. This theory supposes that carrier multiplication in thyristors occurs due to space 


charge effects of holes. [9] 


Figure 8 shows an energy band diagram of a thyristor that is biased slightly below 
the switching point. The p-n hook from Shockley’s theory is formed by the potential bar- 
rier for holes in the P2 layer. The P2 layer collects holes which were emitted from P1 
and carried through N1 by transistor action. In order to understand Shockley’s theory, 
two current components must be considered, electron current coming from N2 and hole 
current coming from Pl. At the switching threshold these currents are beginning to in- 
crease because the potential barriers formed by J1 and J3 are lowering as the junctions 


become forward biased. 


J1 J2 J3 
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Figure 8. | Shockley’s p-n Hook Collector (After [9].). 


13 


The flow of current can be described by following the path of a hole injected from 
P1 into N1 by thermionic emission. This hole would likely cross the relatively short N1 
region without recombining and enter the P2 region by transistor action. The P2 layer 
represents a low potential area for holes, so they will tend to collect there are excess ma- 
jority carriers. The holes then impart a net positive charge to the P2 layer which tends to 
pull in electrons from the N2 layer. Following the reciprocal path, electrons tend to 
gather in the N1 layer, which tends to pull in more holes. This process is thus regenera- 


tive [9]. 


The relative doping level of N2 is much greater than P2; therefore the current 
through J3 is dominated by electrons. Since the width of P2 is about one diffusion 
length, these injected electrons largely transit the P2 layer without recombining. The 
holes injected through J2 into P2 become excess majority carriers. These holes tend to be 
trapped in the P2 region because of the potential barrier formed by J3. This results in 


electron current through J2 that is much larger than the hole current through J3 [9]. 


Clearly, the current multiplication due to the J3 junction is very important to thy- 
ristor operation. However, despite Shockley’s early position that this current multiplica- 
tion is not due to the avalanche process, the idea that avalanche multiplication is involved 
has persisted. Numerous papers and texts dating from the 1950’s-on continue to explain 
current multiplication using the avalanche process [5, 7, 10, 11]. The issue of avalanche 


multiplication is revisited later in this thesis. 


Other models of thyristor operation have been developed, as well. Perhaps the 
most useful of these, at least for appreciating the interactions between the layers, is the 
two-transistor analog. This model is frequently attributed to Moll et al. in 1956 [10]. 
However, the literature shows that in 1952 Ebers mentioned Shockley’s claim to be able 
to get thyristor behavior from two properly connected transistors. Apparently, Shockley 
established the basis for the model but did not put it in writing [12]. The two-transistor 
analog is presented in the following section. 

2. Two-Transistor Analog 

The two-transistor analog considers the thyristor, specifically in this case a 


Shockley diode, as two Bipolar Junction Transistors (BJTs), one npn and one pnp, con- 
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nected as shown in Figure 9. In this figure the two transistors are labeled using both the 
layers of the thyristor (P1, N1 P2, and N2) and the appropriate BJT layer: emitter (E), 
base (B) and collector (C). The subscript pnp or npn is used to distinguish between the 
two transistors. It must be pointed out that this model is useful only in showing where 
the device reaches the switching point under forward bias. It does not explain the opera- 
tion above the switching point. Nevertheless, it is very useful in qualitatively considering 


thyristor switching. 






CATHODE 
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Figure 9. Two-Transistor Analog for the Thyristor. 


When forward bias is gradually applied to the device at equilibrium a small cur- 
rent begins to flow. The emitter layers begin injecting holes from P1 into N1 and elec- 
trons from N2 into P2 by thermionic emission. Because the N1 and P2 regions are ap- 
proximately one diffusion width wide, most of these carriers transit their respective BJT 
bases without recombining and enter in to the collectors. This current tends to forward 
bias both emitter-base junctions (J1 and J3) and tends to reverse bias the base-collector 


junctions (J2). In BJT terms, this biases both transistors in the active mode. 


The holes gathering in the pnp collector (P2) and the electrons in npn collector 
(N1) have no external escape. If these carriers can be seen to move from one transistor’s 
collector to the opposite transistor’s base, they become excess majority carriers in those 
bases. This does not actually involve moving carriers between thyristor layers, only be- 
tween hypothetical transistors. The excess majority carriers in the base can then diffuse 


into the emitter regions: P1 for electrons and N2 for holes. This means that holes dif- 
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fused into N2 will generate electrons injected back into P2, and electrons injected into P1 
will generate electrons injected back into N1. This process is then repeated as the minor- 
ity carriers in the base layers are transmitted into the collectors. Thus the process is re- 
generative [4]. The switching point in the two transistor model is equated to driving both 
BJTs into saturation, or forward biasing all pn junctions. It will be shown later that this 


actually happens at the holding point, not the switching point. 


The current multiplication that occurs at switching threshold can be understood by 
considering the relative doping levels between the emitters and the bases. In both cases 
the emitter is much more heavily doped than the base, or Pl>> NlandN2> P2. Ina 
BJT, emitter efficiency y is defined as the fraction of total current through an emitter- 
base junction which is due to majority carriers in the base. In the case of J1, this would 
be given by [4] 

le ii - e (2.1) 
where /,,, represents hole current through the junction, /,,, represents electron current 


through the junction, and /,, represents the total current through the junction. Another 


expression for emitter efficiency in terms of material parameters is given by [4] 
1 


y= 
14[ De Le No sinh(W/L,) 
D, L; N, ) cosh(W/L,) 


(2.2) 








where Dz, represents diffusivity in the respective layers, Lz represents diffusion lengths 
in the respective layers, Nz represents doping levels in the respective layers, and W 
represents the width of the base layer. If (2.2) is reduced to 

1 


ieee 


E 





ih (2.3) 


by assuming the product of the other terms is unity, then ycan be seen to depend on rela- 
tive doping of the base and emitter. This assumption is made here to demonstrate a point 


regarding current multiplication. However, in most cases the difference between doping 
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levels is much larger than the other differences, so doping level difference is usually the 


primary factor in determining y [4]. 


As an example, where P1 is doped at 1x10'°cm® and N1 is doped at 1x10'°, (2.3) 
gives an emitter efficiency of 0.999. Then rewriting (2.1) to solve for /,,, when /,,, is 


represented by one electron diffusing from N1 into P1, the result is 999 holes injected 
from Pl into N1. This example demonstrated the current gain at J1, and a similar process 


is happening at J3. 


Another thyristor explanation has been created which uses the variation in current 
gain of the two equivalent transistors to define a switching condition. Moll et al. define 


the switching point in terms of the sum of the current gains of the transistors, @,,,, and 
@,,,- The thyristor is in the forward blocking mode when @,,, + @,,,, <1 and above the 
switching point when @,,,+@,,, 21. The current gain for each BJT is not a constant 


value, rather it is a strong function of current density. As transistor current increases, as 
with increasing forward bias, the current gains increase until their sum exceeds unity and 
the thyristor switches [10]. Moll’s explanation is based on a static case. Subsequent 
work has extended this analysis to the dynamic case [13]. Bear in mind that the values 
for current gain in the equivalent BJTs are much lower than typical BJT values. The thy- 
ristor structure does not correspond with ideal BJT structure, so BJT-based parameters 
for a thyristor are bound to be less than ideal [4]. 

3. Junction Capacitance 

The model of thyristor behavior presented thus far has not relied upon fundamen- 
tal semiconductor physics. This is because these models were derived from a qualitative 
point of view. Nevertheless, some aspects of semiconductor physics can be introduced to 
better understand the operation of the thyristor. An important property for thyristor op- 


eration is junction capacitance. 


A pn junction has two capacitance components, depletion capacitance and diffu- 
sion (or storage) capacitance. Depletion capacitance is normally associated with reverse- 


biased junctions since it is the only capacitance component under reverse bias. Deple- 
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tion capacitance also exists in forward-biased junctions, but it is swamped by the other 
component [14]. The source of depletion capacitance can be explained by considering 
the charge remaining in a junction at equilibrium as shown in Figure 10. The depletion 
region consists of two oppositely-charged regions of ionized dopant atoms. Majority free 
carriers, holes on the left and electrons on the right, are then separated from each other by 
the depletion width. This separation of opposite charge resembles a parallel plate capaci- 


tor. 
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Figure 10. Depletion Capacitance. 


As applied voltage V4 across the junction becomes more negative, that is the junc- 
tion becomes more reverse-biased, the distance between the free charges increases and 


the capacitance decreases. Junction capacitance for an abrupt junction is given by [14] 








Guess geENNa 


for V,, <V,. 2.4 
: 2(V;,—Va )(N, +z) ee 


bi 


where A represents cross-section area, g represents elementary charge, € represents the 
permittivity of the material, N, and N; represent ionized acceptor and donor concentra- 


tions on either side of the junction, and V;; represents built-in junction potential, which is 


given by [14] 
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Equation (2.4) is only valid for applied voltages less than the built-in voltage. 


The second capacitance component, known as diffusion or storage capacitance, 
becomes important as the junction becomes forward biased. As the junction applied bias 
V4 becomes more positive, the injection of minority carriers across the junction also in- 


creases. Diffusion capacitance C, is due to the rearrangement of minority carrier den- 


sity [5]. Under forward bias, excess minority carrier concentrations accumulate on both 
sides of the depletion region with a maximum concentration at the boundary. This cre- 
ates accumulations of opposite charge separated by a small distance. This is, of course, 
capacitance. A low-frequency expression for diffusion capacitance at an abrupt junction 
is given by [5] 
oe Aoi Hake Pate ogy Mat (2.6) 
kT 2 2 kT 
where L,,, represent diffusion lengths, kg represents Boltzmann’s constant, T represents 


lattice temperature, and 7, and p,, represent equilibrium minority carrier concentra- 


po 


tions. 


The relationship between depletion and diffusion capacitances is well defined for 
reverse-biased junctions, where only depletion capacitance exists, and strongly forward- 
biased junctions where diffusion capacitance dominates. Less well known is the relation- 
ship between the two components for the low forward-biased case. Casey reports that for 
a junction with built-in voltage of 0.9 V, diffusion capacitance dominates for all applied 
biases greater than 0.8 V [14]. This is shown graphically in Figure 11. Since diffusion 
capacitance is negligible unless the junction voltage approaches the built in potential Vi, 
it can be ignored for low forward-bias and reverse-bias cases. Similarly, junction capaci- 
tance increases at a much slower rate than diffusion capacitance, so it can be ignored for 


strong forward-bias situations. 
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Figure 11. Relation Between Depletion Capacitance and Diffusion 
Capacitance (After [14].). 


This thesis considers only the static case of thyristor operation. The role of junc- 
tion capacitance in static operation is clearly limited but it must play a great role in dy- 
namic operation. The thyristor itself, since it contains three junctions, can be modeled as 
three capacitors in series. Moreover, these capacitances change with bias across the de- 
vice. The interaction of these three junction capacitors, along with the action of the ex- 
ternal capacitor shown in the circuit in Figure 1, are likely to be key to understanding the 


pulsed output from the thyristor detector. 


C, SIMULATIONS USING ATLAS BY SILVACO 

1. Introduction 

ATLAS by Silvaco, Inc. is a powerful, physics-based semiconductor simulator 
program. It predicts the electrical behavior of charge carriers in semiconductor structures 
by finding numerical solutions to user-specified carrier transport equations [15]. 
ATLAS also allows user-selection of a full range of models for physical phenomena such 
as impact ionization, tunneling, etc. The models used in this thesis are explained in the 
following section, along with a brief introduction to how ATLAS was used in the rest of 
this thesis. This is followed by some results obtained from simulations run on a typical 


thyristor. 
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2. Physical Models 
The first step in performing simulations using ATLAS is determining which mod- 
els to use. This section describes the models chosen to simulate the thyristor. These 
models were chosen to simulate relatively small thyristors, with minimum dimensions of 
about | micron, and relatively low switching voltages of about 5V. These models are 
also expected to be valid for larger thyristors with higher switching voltages but that is 
beyond the scope of this thesis. 
a. Drift-Diffusion Transport Model 
The drift-diffusion model of carrier transport explains the motion of 
charge carriers through drift in an electric field and diffusion from high concentrations to 
low concentrations. This model is generally sufficient to explain the operation of devices 
with minimum dimensions greater than one micron [16]. The drift-diffusion model is de- 


rived from the more general Boltzmann Transport Equation (BTE). 


The BTE describes the time-dependent carrier distribution in six- 


dimensional position-momentum space. It is given by [16] 
of. os 
ap te V ETE Vf HCY) (2.7) 


where fis the position-momentum carrier distribution function, 7 and p represent position 
and momentum, F represents the total force acting on carriers, v represents carrier veloc- 
ity and C represents the rate of change of the distribution fdue to collisions. The BTE 
can be derived from conservation of carriers in position-momentum space and from car- 
rier trajectories in position-momentum space [17]. Although the BTE gives an accurate 
representation of carrier distribution, it is not easily solved. In most realistic cases nu- 
merous simplifying assumptions must be made. One set of these assumptions leads to the 


drift-diffusion model. 


The drift-diffusion transport model was derived from the BTE using the 
simplifying assumption that carriers are in equilibrium with the local electric field [17]. 
The simplest version of the drift-diffusion model is based on two equations which de- 
scribe the motion of the two carriers due to electric field drift and concentration diffusion, 


two current continuity equations, and two additional equations describing charge concen- 
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tration and trap density variation. This results in a set of six coupled partial differential 
equations. 

(1) Carrier Motion Due to Drift and Diffusion. The current densi- 
ties due to drift and diffusion add together to give the conduction current density J. The 
two equations, with subscript n and p for electrons and holes, respectively, are given in 


one dimension as [5] 








J, = qltynE + qD, 2 = qu, ne + Not on (2.8) 
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Then the total conduction current is given by Je = Jn + Jp. 

(2) Carrier Continuity Equations. The carrier continuity equations 
are used to account for the flux of free carriers in to and out of an infinitesimal volume. 
Within this volume, the total rate of carrier concentration increase is the algebraic sum of: 
the rate of carriers flowing in to the volume, the rate of carriers flowing out of the vol- 
ume; the rate of carriers being generated within the volume; and the rate of carriers re- 


combining in the volume. The continuity equations are given by [14] 
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where G,,,, and R,,, represent the generation and recombination rates in the material under 
consideration. These equations are valid for consideration of both minority carriers and 
majority carriers. To focus on the minority carriers, which are more important in describ- 
ing carrier transport in bipolar devices, these equations can be modified to represent mi- 
nority carrier continuity. To do this, the current density terms in the continuity equation 
are replaced by the drift-diffusion current density expression. After talking the deriva- 
tives (2.10) and (2.11) reduce to [14] 
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2 
oa = pull, + 1,6 +D, ma +(G,-R,). (2.13) 
In the steady state, defined by d/dt = 0, the electric field term E is 
independent of time. Therefore the partial derivative of E becomes an ordinary deriva- 
tive. The change in the electric field strength with respect to distance can be related to 
the net charge concentration p by Gauss’s Law [14] 


dx E 
where ¢ represents the dielectric constant. The meaning of (2.14) is that when the electric 


(2.14) 


field in a region is constant, that region must also be balanced in charge. Conversely, any 
region not in charge balance must add to or detract from the electric field. Thus in a neu- 
tral region the first and second terms on the right side from (2.12) and (2.13) are elimi- 


nated. The continuity equations in neutral regions simplify to [14] 
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These equations reduce even further in the steady state to 
o 
D,—2+(G,-R,)=0 (2.17) 
ox 
O’P, 
Dis +(G,-R,)=0. (2.18) 


Finally, in the separate case when carrier distribution is uniform, (2.15) and (2.16) reduce 
to [14] 
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(3) Additional Equations. The remaining two equations in the 
drift-diffusion set are described below. The first expresses conservation of carriers in 
trap states. It is given by [18] 


Oi 
5 =(G,-R,)-(G,-R,). (2.21) 
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In most cases, including the device discussed in this thesis, there is no appreciable time 
change in the number of trap states, so (2.21) can be ignored [18]. The last expression 


accounts for the generation of electric field from net charge [19] 


V-E=4[ p-n+Ni-N; |. (2.22) 
E 


This can also be expressed in terms of electrostatic potential V [18] 


-VV = p=4| p-n+Nj-N;, |. (2.23) 
E 
This results in a set of six coupled partial differential equations 
comprised of (2.8), (2.9), (2.10), (2.11), (2.21), and (2.22). This set is usually reduced 
to five by disregarding (2.21). Semiconductor simulator programs using the drift- 
diffusion model simply solve these equations numerically at many points in the device. 
This is the approach used by ATLAS in this thesis [18]. 
b. Impact Ionization 
Impact ionization occurs when free carriers gain sufficient kinetic energy 
from local electric fields to break covalent bonds of bound electrons upon collision. This 
process requires a high electric field over a long enough distance to accelerate the carriers 
to the ionization energy £;. Furthermore, the distance between collisions must also be 
large enough to allow acceleration to the ionizing velocity before a collision removes the 
kinetic energy from the carrier. The general expression for impact ionization rate G is 
given by 
G=a4,J,+a,J, (2.24) 
where @, , represent ionization coefficients. The ionization coefficients represent the 


number of electron-hole pairs generated by a carrier per unit length traveled when the ki- 


netic energy of the carrier is greater than £; [15]. 


ATLAS allows the user to choose from a number of impact ionization 
models. These models can be classified into two main types, local models (which are de- 
pendent on the local electric field) and non-local models (which depend on the kinetic 


energy of the carrier). 
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Local models will normally overestimate the rate of impact ionization 
[15]. This can be explained by considering the motion of a carrier through an electric 
field with a sharply-peaked magnitude above the ionizing threshold. This structure is 
typical of the field at the metallurgical junction of a p-n junction. The local model will 
then calculate the impact ionization across the entire width of the peak. However, since 
carriers require acceleration in that electric field over some relatively large distance be- 
fore they actually exceed ionization energy, the local model will overestimate the ioniza- 
tion rate [15]. Non-local models calculate impact ionization by considering carrier ki- 


netic energy rather than electric field strength. 


A local-type model was chosen for the initial simulations because of the 
overestimation tendency of local models. It will be shown later that impact ionization is 
not strongly involved in small thyristor operation. In certain cases impact ionization does 
not occur at a high enough rate to change the IV curve of the device. The ATLAS model 


selected for this simulation calculates the ionization coefficients using 
a, = Ayexp (-2) (2.25) 


a, = Apexp (-2e) (2.26) 


where Ay p and By p are user-definable parameters and F represents electric field strength 
[15]. 

c. Mobility 

ATLAS also allows numerous mobility models to be used in simulations. 
Simulations of bipolar devices, such as the thyristor, can typically be accomplished using 
two mobility models, concentration-dependent mobility and field-dependent mobility. 
The concentration-dependent model uses a look-up table to assign mobility values based 
on the total doping concentration at a point. These data have been determined through 
empirical observations. Several mathematical models have been derived to describe the 
relationship between concentration-dependent mobility and total doping concentration. 


One useful mode for silicon at 297 K is given by [14] 


2a 


= ol et LR (2.27) 
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The field-dependent model is used to model carrier velocity, including the 





and 


velocity saturation effect at higher electric fields. The field-dependent equations are 


given by [15] 


2, 
Us, (E) = Ly {4 (2.29) 
E -l 
u.(E)=n,,|1+| 2 (2.30) 
p po y 
Sat,p 


where Vsar, np represents the carrier saturation velocity, and //, ,) represents the low-field 


mobility for the material. 

d. Recombination 

Two recombination processes were selected for use in these simulations. 
The first is an impurity concentration-dependent version of the Shockley-Hall-Read 
(SRH) model. The SRH process generates phonons from electron-hole recombinations 
which take place using trap levels which exist in the bandgap of real semiconductors. 
The concentration-dependent version of this model accounts for a decrease in carrier life- 


times in the presence of moderate or higher doping levels. The carrier lifetimes 7, are 


determined by [15] 








T, 0 
= (2.31) 
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where N is the total doping concentration. Low-concentration lifetimes 7,,_,, and pa- 


rameter Nsruyp are user-definable, with ATLAS default values of 0.1 ws and 5x10'° 
cm”, respectively. The low concentration lifetime value is low compared to other au- 
thors. Law compiled the research of several researchers and estimated the low- 


concentration to be 30 ws for electrons and 10 ws for holes [20]. 


The second recombination process occurs through the Auger process. 
This is a non-radiative process that requires three carriers at the same point in order to 
occur. Therefore Auger recombination is likely only in highly doped materials [14]. The 
Auger recombination rate is given by 

R, = AUGN( pn’? — nn? )+ AUGP(np” — pn; ) (2.33) 

where n; represents intrinsic carrier concentration and v and p represent free carrier con- 
centrations. AUGN and AUGP are user definable parameters. The ATLAS default value 
for n, is 1.45x10'° cm” in silicon. This value has been disputed by some authors [4, 21]. 
The results contained herein do not specifically depend on calculations using a precise 
value for n,, therefore the ATLAS default value was used throughout. 

e. Band Gap Narrowing 

Bandgap narrowing is the effect that high doping concentrations can have 
on semiconductor materials. In silicon this is not observed for concentrations below 
about 1x10'* cm* [15]. Slotboom and de Graaf derived an expression for bandgap nar- 


rowing AE’, in silicon BJTs, given by [22] 








2 
AE, = 0.009 nf + (jon +0.5 |. (2.34) 


The effect of bandgap narrowing is the increase in intrinsic carrier concen- 


AE, 
N, = No, {eXp ne (2.35) 
B 


where njo represents the low doping value for intrinsic carrier concentration. 


tration n; given by 
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J ATLAS Basics 
This section discusses the types of simulations used in this thesis and the format 
of the results of those simulations. In doing so, a thyristor structure is introduced to dem- 
onstrate these basics. This structure will also be used in Section D of this chapter which 
employs ATLAS to further investigate thyristor operation. 
a. Input Deck 
The basis for any simulation in ATLAS is a file called the input deck. The 
input deck consists of the following five sections: structure specification, material models 
specification, numerical method specification, solution specification, and results analysis. 
The details of writing the input deck are beyond the scope of this thesis. The reader is re- 
ferred to the ATLAS User’s Manual for detailed description on input deck syntax [15]. 
Appendix A contains the input deck used to create the device in Chapter II. Figure 12 
shows the hypothetical thyristor used through the rest of this chapter. 








Figure 12. Hypothetical Device and ATLAS Input Deck. 


The thyristor used in the initial simulations was a simplified, non-realistic 
structure. The doping levels and layer thicknesses used in this notional device were es- 
tablished by combining example levels provided in the literature [4] which were refined 
through multiple simulation runs to produce a uniformly doped layered device with a 


relatively low switching voltage V,. The non-physical aspect here is the strictly uniform 
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doping profile shown in Figure 13. This ideal step-function is not achievable in actual 


devices, but allows for more clear illustration of the physical conditions here. 


Ji J2 J3 
















Doping 
10cm? 


Doping 
10’ cm? 


Doping 
10'7 cm? 


Doping 
1019 cm? 


| ul 


B 


_4 Lil Lt al 











1 2 3 4 5 6 7 8 
Microns 


Figure 13. Doping Profile and Layer Thicknesses. 


All the simulations performed in this thesis were done in two dimensions. 
This means that many numerical results, particularly those involving currents, must be 
considered as quantity per micron of length in the third dimension. This dimension 
would be into the page. For example, if the current output from a simulation is given as 
100 pA, this really means 100 pA per zm of length into the paper. A 100- zm long de- 
vice with the structural cross-section given produces 100 times the current that the two- 
dimensional model predicts. 

b. ATLAS Outputs 

There are two standard output files which come from a completed simula- 
tion in ATLAS, the log file and the structure file. The log file is used to display one 
quantity against another. The best example of this is the current-vs.-voltage curve, or IV 
curve. Some of the other quantities which can be displayed in a log file display include 
carrier concentrations, carrier mobility, carrier temperature, carrier generation rate and, of 
course, voltage and current. The key to understanding the log file display is that the out- 
put gives the specified quantity at a specified point as another quantity is varied. This 


point can be located anywhere inside or on the surface of the device. 
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In the case of the IV curve, the anode current is measured as the anode 
bias varies. Another example would be measuring carrier generation in the depletion re- 
gion of a photodiode as light intensity varies. The simulated forward-biased IV curve for 
the thyristor described in Figure 13 is shown in Figure 14. The key IV parameters are 
switching voltage Vs of 16.7 V, switching current J; of 1 pA, holding voltage Vz of 0.41 
V, and holding current /;; of 14 pA. These points are found using the conventional defi- 


nition of the point where the tangent is vertical, or dV/dJ =0. Again, the current figures 


are per unit micron in the unspecified third dimension. 
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Figure 14. Example Log File Display of Thyristor [IV Curve. 


The structure file displays a specified quantity in all regions of a device at 
a specified operating point. A good example of this would be showing the electron dis- 
tribution at thermal equilibrium. Other quantities that can be shown in a structure file 
display include electric field, current density, band structure, and many others. The pri- 
mary structure file display is a multi-color plot where variations are shown by changes in 
colors. The structure file display also includes a feature called the cutline. This feature 
allows a more accurate numerical-scale display of the desired quantity along a specified 


line through the device. An additional feature of the cutline display is that multiple quan- 
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tities can be displayed on the same plot. These features are demonstrated in Figure 15 
which shows a basic structure file display along with the cutline display. The cutline, 
which is taken along a vertical line through the device and plotted on a horizontal line, 
shows the donor concentration and the electron concentration on a logarithmic scale. The 
plot shows that the logarithm of the electron concentration in p-type regions approxi- 


mately agrees with that predicted by the pn product [14] 
i= (2.36) 


Figure 15 also shows that electrons move away from high concentrations at the boundary 


with a low concentration, as is expected via diffusion. 
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Figure 15. |. ATLAS Structure Display and Corresponding Cutline Display. 


c. Simulation Overview 

ATLAS allows the user to choose from a range of simulation types, such 
as DC, AC, or transient, and to perform them under a diverse range of operating condi- 
tions such as bias, current, temperature, and incident light. The simulations contained in 
this thesis are generally DC analyses conducted both with incident light and without. In 


each case the temperature is assumed to be 300 K. 
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Light in ATLAS can be modeled either as a single wavelength source or a 
multi-spectral source. The model used in this thesis has the same shape as the Air Mass 
Zero (AMO) solar spectrum [23]. The intensity of the spectral components varies with 
the overall intensity, which is specified in W/cm*. This light model was used for conven- 
ience. Future work should include characterization using a intensity distribution which is 
uniform over the desired wavelength range. The relative intensities of the spectral com- 


ponents used for simulations in this thesis are shown plotted in Figure 16. 
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Figure 16. Relative Intensity of Spectral Components. 


This section has introduced the ATLAS simulation program and how it 
will be used in the remainder of this thesis. The models of physical processes assumed to 
be at work in the thyristor were presented along with how those models are presented in 
ATLAS. The following section uses ATLAS to explore thyristor operation in detail and 


to compare the results with operational theory. 
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D. SIMULATION RESULTS 

This section analyzes the results of ATLAS simulations conducted on the device 
shown in Figure 12. These results, along with analysis, are used to observe and describe 
changes in the physical state of the device at different operating points. 

1. Switching Conditions 

This section briefly shows the results of the simulations on the example thyristor 
in the vicinity of the switching point. The goal here was to establish the physical condi- 
tions which occur at this point. Figure 17 shows the change in carrier concentration when 


biased at the switching threshold. 
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Figure 17. Carrier Concentrations (em”®) at Switching Threshold. 


Qualitatively, Figure 17 describes exactly what is expected at the junctions. J1 
and J3 are forward biased, as can be seen in the accumulation of minority carriers at the 
forward-biased junction. For example, the increase in electrons in P1 and the increase in 
holes in N2, both near J1, are characteristics of a forward-biased junction. Conversely, 
the decrease in majority carriers near J2 shows that this junction is reverse biased. This 


depletion of carriers occurs as the electric field developed at J2 spreads further into N1 


a3 


and P2. Notice that the carriers are depleted well into the N1 region at switching thresh- 


old. 


The switching threshold represents the maximum bias across the thyristor, and J1 
and J3 are becoming more forward biased as current increases. This means that the junc- 
tion potential at J2 should also be a maximum at the switching point. This is shown in 
Figure 18 which plots the current against the voltage across J2. Note that since the volt- 
age is “double valued,” the data was generated in ATLAS by increasing the current over 
the desired range and computing the resulting voltage. This plot format of current vs. 


voltage is used throughout this thesis. 
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Figure 18. Junction Voltage Across J2 as Current Increases. 


The maximum reverse bias across J2 also means that the electric field across J2 
will be amaximum. The electric field between J1 and J2 is shown in Figure 19. This 
shows that the electric field reaches a peak strength of 2.1x10° V/cm at J2. The electric 


field extends farther into N1 because of the lower doping level than in P2. 
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Figure 19. Electric Field in the N1 Layer. 
Figure 20 shows the voltages across the junctions J1 and J3. These plots confirm 
that the junctions are becoming more forward biased as current through the thyristor in- 


creases. 


























Figure 20. Junction Voltage Across J1 and J3 as Current Increases. 


2. Impact Ionization 
This section revisits the role of impact ionization in thyristor operation. Recall 


that Shockley asserted very early on that the current multiplication in a thyristor was not 
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due to the avalanche process. His reasoning was that the current multiplication was ob- 
served at low enough bias voltages to make the avalanche rate extremely low. Despite 
Shockley’s early position that this current multiplication is not due to the impact process, 
the idea that avalanche multiplication is involved has not diminished. The simulation re- 
sults in this section will show that while impact ionization may occur in some thyristors, 


it does not occur in all devices. 


First, ATLAS simulations are used to show that thyristors can operate without 
impact ionization. The simulated structure used here uses the same doping profile as 


used previously, but the thicknesses of N1 and P2 are set to 1 um. The doping profile is 


shown in Figure 21. 
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Figure 21. Impact Ionization Demonstration Doping Structure. 


ATLAS was then used to generate the simulated IV curve for this device, once us- 
ing the default impact ionization parameters and once setting those parameters to zero. 
The later case is equivalent to turning of the impact ionization model. Figure 22 shows 
the two IV curves. The curve on the right was generated using the impact ionization 
model with the default parameters used for (2.25) and (2.26), which are described below. 
The curve on the left was generated with these parameters set to zero. The curves are 


identical. 
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Figure 22. IV Curve Comparison, With and Without Impact Ionization Modeling. 


In order to investigate this a bit further, the maximum electric field strength for 
this device was found to be 1.1x10° V/cm at J2, as shown in Figure 23. 
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Figure 23. Electric Field Through the N1 Layer. 


The mathematical models for impact ionization used in these simulations, which 


are repeated here from (2.25) and (2.26), are given by 
B 
a, = A,exp (-2: G37) 
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a, = Ayero(-4E (2.38) 


where EF represents the local electric field magnitude. The default parameter values are 
Ay = 7.03x10° cm", By = 1.23x10° V/em, Ap = 7.03x10° cm", and Bp = 1.69x10° V/cm. 
The parameter Ay p can be interpreted as the maximum number of carriers that can be 


generated per unit distance in a very strong electric field (1.e., E >> By ») while By p 


represents a scaling factor related to where the electric field becomes very strong [24]. 
Equations (2.37) and (2.38) can then be interpreted as meaning that, unless the electric 
field E approaches the values given by By or Bp, the exponential term will dominate and 
drive the ionization rate very low. Comparing the maximum electric field in this device 
to these parameters indicates that this model will show a very low rate of impact ioniza- 


tion. 


The best values to use for the parameters Ay.p and By p are uncertain. However, 
Maes et al. [24] conducted a review of the literature related to impact ionization in sili- 
con. Their results, while inconclusive as far as actual parameter value valid for low field 
strengths, indicate that the rate of impact ionization would be very low in this device. 
Using the simulated peak field strength of 1.1x10° V/cm, Figures 24 and 25 indicate the 
values for the ionization rates a, and @, would be 2 cm" and less than 1 cm”, respec- 
tively. Again, this indicates that impact ionization is not a factor in the operation of this 


thryristor. 
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The previous example shows that a thyristor can be simulated to operate properly 
without impact ionization. However, this does not mean that impact ionization does not 
occur, just that it is not the primary current multiplication mechanism at work. As an ex- 
ample of a device where impact ionization does have some impact on the IV curve, con- 


sider the device introduced in Figure 12 at the beginning of this chapter. The maximum 
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electric field strength from Figure 19 is 2.1x10° V/cm, which is still well below the 
strong field threshold set by parameters Ay and Ap of 7.03x10°cm". As Figure 26 
shows, the IV curve does show a small effect between simulations with the impact model 


enabled and disabled. 
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Figure 26. IV Curve Showing Impact Ionization Effects. 


The result of this section is that impact ionization is not the primary cause of the 
current multiplication in a thyristor. While impact ionization may occur in a thyristor, it 
is not critical to the operation of the device. 

3. Holding Conditions 

Recall that the two-transistor model for thyristor operation described the switch- 
ing point by finding where both BJT’s become saturated. It was previously shown that 
this does not happen at switching but can be used as criteria for establishing the holding 
point. At some anode current value, J2 must become less reverse biased than it is at equi- 


librium. This can be carried out by ATLAS using two approaches. 


The first approach is a graphical analysis of structure plots of the electric field 
near J2. As J2 becomes slightly forward biased, the electric field will be of a slightly 


lower magnitude and of slightly smaller lateral extent into the N1 layer. Figure 27 shows 
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the electric field near J2 simulated at equilibrium and at the holding point as defined in 
Figure 14. This holding point definition was determined by finding the point on the IV 


curve with a vertical tangent; in Figure 14 it was found to be 14 pA. 
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Figure 27. Electric Field Near J2 at the Holding Point. 


Figure 27 clearly shows that the electric field at the holding point, shown in blue, 
is significantly less than at equilibrium, shown in red. This indicates that J2 becomes for- 
ward biased at a lower anode current than that defined by finding the vertical tangent to 


the IV curve. 


Further simulations were performed until a more accurate value could be deter- 
mined to define when J2 becomes forward biased. Figure 28 shows the electric field near 
J2 again, but this time the blue line represents an anode current of 5.7 pA. In this case, 
the electric field has a slightly lower peak magnitude and extends slightly less distance 
into N2 than at equilibrium. This is a better definition for the holding point than the ver- 


tical tangent to the IV curve. 
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Figure 28. Electric Field Near J2 at the Holding Point. 


The second method of finding where J2 becomes forward biased is to find the an- 
ode current where the junction potential across J2 (Vj) is equal to the built-in potential 
Vi across the junction. The built-in potential of a junction was given by (2.5), which is 


repeated here for convenience 





b 
Nn. 


V7= eae in( Be (2.39) 
q i 
Equation (2.39) applied to J2 gives a built-in voltage of 0.756 V. Figure 29 shows a dif- 
ferent view of the data contained in Figure 18 with this built-in potential plotted to de- 
termine the corresponding anode current. This value, indicated by the dot, is approxi- 


mately the same as determined using the prior method. 
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Figure 29. Anode Current vs. J2 Junction Potential. 


Figure 30 shows the location of the holding current as determined above com- 
pared to the value determined by the vertical tangent. The graphical solution occurs at a 
much lower anode current, and therefore does not occur at the minimum applied voltage 


across the entire thyristor. 
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Figure 30. Close Up View of IV Curve. 


4. Variation of parameters 
This section discusses a series of simulations performed on the basic device struc- 


ture already introduced, but with one of the key parameters varied while the others are 
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held constant. This will allow observation of the effect that parameter changes have on 
the IV curve. The device parameters that will be varied are the doping concentration of 
the two BJT emitters, Pl and N2, and thickness of the two BJT bases, N1 and P2. In 
each plot the IV curve for the device as already described will be given as a reference. It 
will be indicated in the legend of each figure by the word reference attached to the file 
name. 

a. Variation of Emitter Doping 

The P1 layer is described in the two-transistor model as the emitter for the 
pnp transistor and the N2 layer is the emitter for the npn transistor. Equation (2.3) de- 


fined the emitter efficiency yof a BJT with the assumption that variations in ywere de- 


termined primarily by changes in the relative doping of the base and emitter. This section 
continues that assumption to explore the changes to the simulated IV curve as the doping 
levels are varied in the two BJT emitters. The doping levels were changed separately to 
allow observation of the effects on the switching voltages and holding currents independ- 


ently. 


Figure 31 shows the change in the IV curve as P1 layer doping is in- 
creased. At the low end of this range, emitter efficiency at J1 is very low, and approaches 
unity at the high end of the range. The legend in Figure 31 shows the name of the 
ATLAS output file which corresponds to the various doping concentrations in the P1 
layer, which are given in scientific notation form. Thus, le15, shown in red represents 
the IV curve simulation result when P1 is doped to 1x10" cm”. This notation is used 


through this section when doping levels are varied. 
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Figure 31. Anode Current vs. Bias as the P1 Doping Parameter is Varied. 


Figure 31 shows that high emitter efficiency at J1 enhances the completion 
of the turning-on process by lowering the holding current. The lowest holding current 
occurs at the highest emitter efficiency. Decreasing the doping level decreases the emit- 
ter efficiency, which also increases the holding current. The lower doping also increases 


the switching voltage. 


Figure 32 shows the same doping changes applied to the N2 layer. As 
with the previous case, lower N2 doping results in lower emitter efficiency at J3, higher 
holding current, and higher switching voltages. However, unlike the previous case, the 
lowest holding current does not occur at the highest emitter efficiency. This suggests 
that, in the case of N2 doping, there exists an optimal doping level to enhance thyristor 
turn-on. This optimum doping is not equal to the maximum feasible doping as in the 


previous case. 
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Figure 32. Anode Current vs. Bias as the N2 Doping Parameter is Varied. 


As shown in Figures 31 and 32, the lowest emitter doping cases, shown in 
red, resemble the reverse-bias breakdown curve of a diode, which is caused by impact 
ionization. The higher bias at the switching point for the low emitter doping cases also 
means that impact ionization happens at a higher rate. Figure 33 shows an additional 
simulation performed using the N2 emitter doping level of 1x10'°cm”. In this figure, the 
red line represents the simulation with the impact model enabled and the blue line repre- 
sents the simulation with the impact model disabled. This indicates, again, that impact 


ionization is not the dominant factor in thyristor operation. 
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Figure 33. TwoIV Curves With and Without the Impact Ionization 
Model Enabled. 


The increased switching voltage for low emitter doping can be explained 


using the theory that switching happens when the sum of the current gains @,,,, @,,,, are 


equal to or greater than one. Low emitter doping will reduce the current gain, thereby in- 
creasing the switching current. In the lowest doping case for N2, the switching current is 
an order of magnitude higher than the higher doped cases. This indicates that low emitter 
efficiency is the primary explanation for the higher switching voltage and current for low 


emitter dopings. 


From the consideration of designing a thyristor detector, Figures 31 and 32 
suggest that lower emitter efficiency is desired to enhance the operation of a thyristor- 
type detector. This will increase the holding point to make it less likely that the turn-on 
process is completed. However, the increase in switching voltage may make the holding 
current improvement unacceptable. 

b. Variation of Base Thickness and Doping 

In this section the thicknesses and doping levels of the N1 and P2 layers, 
which are the bases of the pnp and npn transistors in the two-transistor model, respec- 
tively, will be varied. The IV curve as the P2 base thickness is varied is shown in Figure 


33. Increases in P2 thickness have the predictable effect of increasing the switching 
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voltage and holding current. The basic shape of the IV curve, however, remains the same 


as P2 thickness varies. The legend for Figure 34 shows the names of the ATLAS output 


files which correspond with the simulation results for various layer thicknesses, which 


are specified in microns. Thus, 1_5.log corresponds to the simulated IV curve when the 


thickness of the P2 layer is 1.5 4m. This convention is used throughout this section the 


for layer thickness variation plots. 
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Anode Current vs. Bias as the P2 Thickness Parameter is Varied. 


Figure 35 shows the changes to the IV curve as the P2 doping level is var- 


ied. The effect of increasing the doping level of P2 is similar to decreasing the doping 


level of N2; it reduces the emitter efficiency. Figure 35 shows, as does Figure 33, that 


lower emitter efficiency results in higher switching voltage and holding current. The in- 


crease in switching voltage with increased thickness can be explained by considering that 


fewer electrons injected through J3 will successfully transit the P2 when it is thicker. In 


BJT terms, a thicker P2 layer results in a lower base transport factor [4]. 
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Figure 35. 


Next, the same variations were performed on the N1 doping and thickness. 
Figure 36 shows the simulated IV curve as N1 thickness is varied. It is immediately evi- 
dent that the operation of the thyristor is much more sensitive to changes in N1 than in 


P2. Notice that for IV curve associated with N1 thickness of 0.5 “zm the device does not 
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Figure 36. Anode Current vs. Bias as the N1 Thickness Parameter is Varied. 
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Figure 37 shows the change in the simulated IV curve as N1 doping is var- 
ied. Again, the switching voltage is strongly influenced by the changes in the doping 
level of N1. 
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Figure 37. Anode Current vs. Bias as the NI Doping Parameter is Varied 


Figures 36 and 37 show that the thyristor switching point is primarily de- 
pendent on the doping and thickness of the N1 layer, with a secondary dependence on the 
P2 doping. These parameters can be used to determine the switching point for a thyristor, 
and the remaining parameters, including the doping of P1, P2 and N2, can be used to 
shape the IV curve as desired. 

=) Optimal Thyristor Detector Parameters 


The results of the previous section can be used to develop a few conclusions re- 
lated to optimizing detector operation of a thyristor. These results are necessarily general 


relationships between layer parameters. 


The first goal of a thyristor designed to produce a pulse-mode output is that it 
must not complete the turning-on process. That is, the anode current must not exceed the 
holding current. This can be achieved most effectively by lowering the efficiency of the 
emitters in the two-transistor model, Pl and N2. This can be done by decreasing the dop- 


ing of the emitters, P1 and N2, or increasing the doping concentration of the P2 base. 
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Lesser effects can be achieved by increasing the thickness of the N1 or P2 layers, or in- 


creasing the doping concentration of the N1 layer. 


The switching voltage can also be controlled by varying the parameters. The pri- 
mary determinant in the switching voltage is the doping and thickness of the N1 layer. 
These parameters can be used to select a switching voltage slightly above the intended 


biasing voltage for the thyristor. 


E. CHAPTER CONCLUSION 

This chapter has introduced the thyristor device and described the static theory of 
thyristor operation. ATLAS simulations were introduced and applied to demonstrate thy- 
ristor operation in different operating modes. It was shown that the switching point volt- 
age is determined primarily by the thickness and doping of the N1 layer and the doping of 
the P2 layer. The final section of the chapter showed how device parameters can be used 
to modify the shape of the IV curve. The following chapter will apply some of these 


items in the effort to design a thyristor to be used as an optical detector. 
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Hil. DEVICE DESIGN AND SIMULATION 


This chapter discusses the design and simulation of the photodetecting thyristor. 
As originally proposed, this thesis would have stopped at the computer simulation stage. 
However, the author was awarded a research Fellowship from the Space and Naval War- 
fare Systems Center, San Diego, to extend the research effort to include fabrication of a 
device to test the concept. Because the late funding of this project, the design has been 
submitted for fabrication but the fabrication is not yet complete. This chapter presents an 
overview of the process used to create the device design, an introduction to the design 
created, and the results of simulations performed on the device. The device designed is 
not an ideal detector structure; rather it is a functional design which can be built using an 


existing fabrication process. 


There are two procedures used to create the device design. The first is creating 
the device in ATLAS to optimize placement of the thyristor layers. This is done using 
the constraints of the fabrication process parameters and the layout design rules. The 
second procedure converts the design to a layout file which is used to specify how the 
device is to be fabricated. The second step was done using L-Edit by Tanner, Inc. Be- 
fore the design could be started, it was necessary to select a fabrication process. The 


steps taken to select the process and design the device are described below. 


A. PROCESS SELECTION AND PARAMETER ESTIMATION 

There are very limited options available to fabricate small runs of integrated cir- 
cuit devices at an affordable cost. Perhaps the most widely used fabrication service for 
small lots is provided by the MOSIS service. MOSIS is not a manufacturer; rather it is an 
organization that compiles multiple small-run projects submitted by companies, govern- 
ments, and universities onto a single wafer which is then fabricated by a semiconductor 
fabrication company. MOSIS works with several different manufacturers allowing for a 
wide variety of available processes. The result is that MOSIS provides a method for get- 


ting small-lot projects manufactured using a specific process at a reasonable cost. 
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The considerations used in choosing a fabrication process include minimum fea- 
ture size, cost, and the number and types of layers available to the circuit designer. This 
project does not involve particularly small features or excessive numbers of layers, but it 
does require that the process permit bipolar devices. The bipolar requirement proved to 
be the real limit in selecting the process. MOSIS offers only one such process, the 1.50- 
uum ABN CMOS process (ABN) by AMI Semiconductor, Inc. As with any commercial 
process, the circuit designer has no control over the doping concentration and thickness 
of the various layers. This required estimation of the process parameters to permit simu- 
lations. 

1. The AMI ABN Process 

The ABN process allows up to two metal layers, two polysilicon layers, an n well 
and a p base, all on a p substrate. This process was designed to permit a vertical npn BJT 


as shown in Figure 38. 


Vertical Transistor 
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nwell (Collector) 
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Figure 38. Cross-Section View of the ABN Process Vertical 
BJT Structure (After [25].). 


This basic scheme was modified slightly to create the pnpn thyristor needed for 
this project by replacing the p+ collector contact with an n+ layer. This allowed the lay- 
ers of the thyristor to be arranged as shown in Figure 39. Rather than a vertical BJT, the 


layout becomes a lateral thyristor. 
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Figure 39. Cross-Section View of the Thyristor Structure 
using the ABN Process. 


2. Parameter Estimation 

One of the difficulties encountered using a commercial fabrication process is that 
the doping levels, doping profiles, and layer thicknesses are not available to the designer 
because they are considered proprietary data. The designer is then tasked to design a de- 
vice based on fixed, unknown parameters. This is not usually a problem for most circuits 
since they typically use manufacturer-provided circuit elements. These are sub-circuits 
and devices which are known to be designed correctly to perform a specified operation. 
Then the circuit designer’s job becomes connecting these known-good elements into a 


larger circuit to accomplish a task. 


This project, however, is not based on these predetermined device layouts. 
Rather, the goal here is to create an entirely new device whose operational characteristics 
depend on the unknown parameters. Specifically, this includes the doping concentrations 
of all the layers and the thicknesses of the N1 and P2 regions. The importance of these 


parameters meant that estimates of the critical parameters were needed. 


Doping concentration in the various layers can be derived by first determining the 
resistivity 9 of each region. This can be done by using 
PH=R;x, (3.1) 
where Rs represents the sheet resistance of the material and x represents the thickness of 


the material. Sheet resistance is non-proprietary data which can be obtained from MOSIS 


aD 


for previous manufacturing runs [25]. These values are expected to remain nearly con- 
stant for future fabrication runs. The value of Rs for the thyristor layers, averaged over 
the last twenty fabrication runs, is given in Table 1. To complete the calculation of resis- 


tivity 9, the layer thickness x must be known or estimated. Since these thicknesses are 


not known, it was necessary to make estimates based on comparison with similar CMOS 
processes [26, 27]. Once the resistivities of the layers were estimated, the doping con- 
centrations were estimated graphically using Figure 40. Estimated layer thickness, resi- 
tivity, and doping concentration are given for each of the layers used in the thyristor in 


Table 1. 
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Figure 40. _ Resistivity vs. Impurity Concentration (From [14].). 
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Region | Average Rs Estimated Resistivity | Doping 
(Q/square) | Thickness (wm) | (Q-cm) (cm*) 
































Pl 76.5 0.3 0.0022 4x10!” 
NI 1598 4 0.64 7x10" 
p2 2265 2 0.45 2x«x10'° 
N2 52.5 0.3 0.0016 4x10” 





Table 1. Critical Parameters of Device Layers. 


B. DEVICE DESIGN 

Before the ATLAS model could be created it was necessary to understand the 
manufacturer’s restrictions and limitations on dimensions and placement of the various 
layers for the process chosen. These restrictions are specified as design rules. This sec- 
tion describes how the basic device was created within the restrictions imposed by the de- 
sign rule set. 

1. Design Rule Non-Disclosure Agreement 

The actual design rules for the ABN process are protected by a non-disclosure 
agreement with the manufacturer. Therefore the actual dimensions of the device as sub- 
mitted for fabrication will not be given in this thesis. The simulations and results de- 
scribed in this chapter were obtained using parameters that are not derived from the pro- 
tected data. Rather, the design is presented using the design rules for the scalable CMOS 
(SCMOS) process, which are available on the MOSIS website [25]. 

a SCMOS Design Rules 

This section introduces the design rules used to create the simulation of the thyris- 
tor device. As discussed above, these are not the design rules used for the actual fabrica- 
tion of the device because of the non-disclosure agreement, rather they are the SCMOS 
rules. Only a small portion of the entire SCMOS rule set will be introduced. Specifi- 
cally, only those rules needed to describe the device for simulation using ATLAS are dis- 
cussed. This excludes the rules for contacts, metal layers, interconnects, and bonding 
pads since these are not used for the simulation. The full set of SCMOS design rules is 


available on the MOSIS web site [25]. 


The primary feature of this design is a thyristor structure created using the basic 


structure of the p base layer inside the n well layer. Figure 41 shows a graphical repre- 


af 


sentation of the design rules pertaining to these layers. The numbers in Figure 40 are ref- 
erences to the rules which specify minimum dimensions allowable in the SCMOS rule 
set, which are listed in Table 2. The units in Table 2 are given in terms of the unit 
lambda. Lambda is a parameter used in SCMOS to allow uniform scaling of all device 
dimensions by the same factor. This allows the same design layout to be used in different 
technologies without resizing the layout. The layout for this project was not drawn using 
lambda-based rules but some value is needed to convert the SCMOS rules to a device 
which can be simulated. For the purposes of this chapter, lambda will be equal to 


0.8 um. The design rule set will be used as shown, except for the N+ Select region near 


the top. As explained earlier, this will be replaced by a P+ Select layer to obtain the pnpn 


thyrsitor structure. 


N+_Select 





Figure 41. _SCMOS Design Rules Diagram for NV Well and P Base 
Layers (Top View) (From [25].). 
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Rule Description Lambda 
16.1 Contact 2x2 
16.2 Minimum N+_Select overlap of contact 3 











16.3 Minimum PBase overlap of N+_Select 





16.4 Minimum spacing between N+_Select and P+_Select 





16.5 Minimum PBase overlap of P+_Select 





2 

4 

2 

16.6 Minimum P+_Select overlap of Contact 2 
16.7 Minimum N_Well overlap of PBase 6 
4 

2 

3 

2 








16.8 Minimum spacing between PBase and Active 





16.9 Minimum Active overlap of Contact 
16.10 | Minimum N_Well overlap of Active 


16.11 | Minimum P+_Select overlap of Active 
Table 2. | SCMOS Design Rules for N Well and P Base Layers (After [25].). 























nF ATLAS Input Deck Creation 

Once the design rules are known, they were then used to create the cross-section 
of the device to be used in ATLAS simulations. This section discusses the layers used in 
the ABN process, which are shown in Figure 41, and how they are used to create the thy- 


ristor device shown in Figure 42. 
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Figure 42. Thyristor Cross-Section View. 
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a. Substrate 

The substrate used in this device is doped lightly p type silicon. It is not 
used in the ATLAS simulations since it is assumed to provide isolation between adjacent 
devices on the actual chip. The substrate is not shown in either Figure 41 or 42. 

b. N Well 

The N_Well contains the entire thyristor device structure. It also acts as 
the N1 layer in the thyristor device. As discussed earlier, the doping level is estimated to 
be 7x10" cm”, and the layer thickness is estimated to be 4.0 um. The width of the n 
well, shown in red in Figure 42, is 50 4m, while the length is either 110 “mor 510 wm. 

c. PBase 

The P Base layer, shown in dark blue in Figure 42, serves as the P2 layer 
for the thyristor. The doping level is estimated to be2x10'° cm” and the layer thickness 
is estimated to be 2.0 um. The P Base layer serves as the P2 layer for the thyristor. It 
also contains the contact to the P2 layer, shown in Figure 41 as P+ Select, and the N2 


layer, shown as N+ Select. The width of the P Base layer 22.8 um. Design rule 16.7 
specifies that the P Base must be at least 6 lambda (4.8 zm) from the edge of the 


N_Well edge. 

d. N+ Select 

The N+ Select layer, shown in light blue in Figure 42, is used only for the 
N2 layer in the thyristor. The doping level of this layer is estimated to be 4x10" cm”, 
and the layer thickness is estimated to be 0.3 4m. The cathode is attached to this layer by 
a Contact layer. Design rule 16.1 states that the contact must be 2 lambda square and rule 
16.2 states that the contact must be 3 lambda from the edge of the N+ Select layer. The 


same contact size is used in all regions in this process. 


Because light is detected by being absorbed in the junction formed by the 
N Well and the P Base layers, it is important that the contacts and the overlying metal 
layers be placed as far away from the junction as possible. This was accomplished by 
making the N+ Select layer wider than the minimum width and moving the contact away 


from the junction. This technique is also applied to the cathode contact to the P1 layer. 
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e. P+ Select and Active 

The P+ Select layer is used to provide the gate contact to the P2 thyristor 
layer, which is formed by the P Base layer, and to create the Pl layer. The doping level 
of this layer is estimated to be 4x10"? cm” and the layer thickness is estimated to be 


0.3 4m. 


The P+ Select region used for the gate contact is located inside the P Base 
layer. It is shown in orange in Figure 42. Design rule 16.4 states that the P Base must be 
placed at least 4 lambda from the N+ Select layer and rule 16.6 states that it must be 
placed at least 2 lambda from the edge of the P Base. The standard contact size is used 


and it must be placed at least 2 lambda from the edge of the P+ Select. 


The P+ Select region used to create the P1 region requires a more careful 
description. The P+ region, shown in green in Figure 42, actually corresponds with the 
Active layer in the ABN process. As shown in Figure 41, the Active layer lies within the 
P+ Select layer. The ion implantation process, which is the method of doping the region, 
only takes place in the active area. In other words, the portion of the P+ Select region 
outside of the Active region retains its original doping. In this case, it retains the doping 
of the N Well region and only the Active region is doped at the P+ level of 4x10" cm”. 
This can be seen in Figure 42 as the gap between the P+ Select and the P Base. This gap 
then defines the width of the N1 layer, which is the primary factor in determining the 
switching point of the thyristor. In the case of the SCMOS design rule, the minimum N1 


width is 4 lambda, or 3.2 4m. This represents the basic device design. 


Because of the relatively large chip area available and the small size of the 
device, three other variations were also designed, both for comparative simulations and 
for fabrication. This allows multiple devices to be built on a single chip. The N1 widths 
of the variant designs are 4.2 um, 5.2 4m, and 2.2 44m. The last N1 width listed violates 
the design rules for the process. This may mean that the fabricated device may not per- 


form as expected, but this will not affect the operation of the other devices on the chip. 


The design rules which apply to the Active region include the following; 
16.8, which specifies a minimum separation from the P Base of 4 lambda; 16.9, which 
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specifies the minimum spacing of the contact from the Active edge of 2 lambda; 16.10, 
provides minimum spacing of the Active from the edge of the N Well of 3 lambda; and 
16.11, which provides minimum spacing of the Active from the edge of the P+ Select of 
2 lambda. 

4. ATLAS Simulation Results 

The design rules discussed in the previous section were used to create an ATLAS 
input deck to simulate the IV curve. The structure for this device was shown in Figure 
42. Figure 43 shows the IV curve for the four variations of the basic device, with the N1 
width indicated in microns in the legend. As expected, the structures with the wider N1 


width have a higher switching voltage. 
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Figure 43. Simulated IV Curve of Fabricated Device. 


Figure 44 shows a structure plot displaying the simulated electric field in the 


3.2 um device biased at the switching threshold. The peak electric field magnitude at J2 


is 65,000 V/cm. According to Figures 24 and 25, impact ionization occurs at a very low 


rate. 
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Figure 44. Electric Field in N1 When Biased at the Switching Threshold. 


C. DEVICE LAYOUT 

Once the design was verified using ATLAS simulations, it had to be converted to 
a useable layout for fabrication. This was done using L-Edit by Tanner, Inc. Figure 45 
shows the top view of the basic device layout. The cross section view shown in Figure 
42 could be taken along a vertical line through Figure 45. The ABN process layers are 
pointed out in this figure, along with the thyristor layer names. The length of the basic 


device, which would be along a horizontal line in Figure 45, was chosen to be 110 4m. 





N+_(N2) 







P Base (P2 


P+ (N1) 
Active (P1) 














Figure 45. Top View of Device Layout. 


The surface area of the chip produced by the ABN process is 2.2 mm square. 


This, of course, is much larger than the area of the device in Figure 45. This allowed 
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multiple devices with different lengths to be placed on the chip. The chip design submit- 
ted for fabrication consists of sixteen separate devices. Each of the four N1 width variants 
are designed with two lengths, 110 “mand 510 wm, for a total of eight distinct device de- 
signs. Two copies of each of these eight designs are included on the chip. A top-view of 


the final chip design is shown in Figure 46. 








Figure 46. Layout of Chip as Submitted for Fabrication. 


D. CHAPTER CONCLUSION 
This chapter introduced the design of the device created and the methods used to 


create that design. ATLAS simulations were used to validate the operation of the device 
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using the principles introduced in Chapter II. The following chapter briefly summarizes 


the work performed for this thesis and makes recommendations for further research. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


This thesis has conducted an in-depth investigation into static thyristor operation. 
The emphasis of this work was to develop a better understanding of thyristor operation to 
enable further work toward designing a thyristor structure for use as a pulse-mode optical 
detector. A key result is the identification of desirable doping concentration relationships 
between the layers needed to lower the emitter efficiency, thereby making the thyristor 


more resistant to completely turning on. 


This thesis is the first of several research efforts needed to advance the idea of a 
thyristor-based detector. The next logical step will be to test and analyze the thyristor 
device described in Chapter III. The previous experiments using a thyristor as a detector 
[2], described in Chapter I, were performed on a device with many unknowns. This in- 
cludes not only the layer structure and doping concentrations of the thyristor but also how 
the device actually sits in the package. Those experiments were necessarily limited to 
qualitative results. Although the new device is not an ideal detector structure, and some 
of the parameters can only be estimated, meaningful quantitative results which can be 


used to improve the structure are expected. 


Another area of research is extending the analysis of thyristor operation into the 
dynamic regime. Of particular interest are the physical changes in the device which pro- 
duce the pulsed output. It is expected that the change in junction capacitance as current 
increases is strongly involved in the pulse generation. ATLAS has the capability to per- 


form transient simulations, so should continue to be of use in this effort. 


Finally, work needs to be done to develop an optimal detector structure. The de- 
vice designed in this thesis has many shortcomings as a detector because of the extremely 
limiting design rules for the process used. Other fabrication processes may have less re- 


strictive rules which will permit a larger area for useful optical absorption. 
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APPENDIX A. ATLAS INPUT DECKS 


The following input deck code was used to generate the simulated IV curve for 


the basic device as described in Chapter III. 


GO ATLAS 


# Sets variable used to input the doping concentrations for the thyristor lay- 
ers 

SET Pl DOPING=4.0e19 

SET N1_ DOPING=7.0e15 

SET P2 DOPING=2.0e16 

SET N2_DOPING=4.0e19 


# Defines variables to be used to define to corners of the various layers 
SET X01=07.2 

SET X02=09. 
SET X03=11. 
SET X04=12. 
SET X05=16. 
SET X06=20. 
SET X07=20. 
SET X08=22. 
SET X09=26. 
SET X10=28. 
SET X11=30. 
SET X12=34. 
SET X13=37. 
SET X14=38. 
SET X15=40. 
SET X16=42. 
SET X17=50. 
SET Y01=0.3 
SET YO2=2 
SET Y0O3=4 


OMWFOANWADWDFAWADDODOAAHNN 


# Defines the mesh used carry out the numerical solution by ATLAS. 
MESH 

-MESH LOC=0.0 SPACING=1 

.MESH LOC=$X01 SPACING=0.5 

.MESH LOC=$X02 SPACING=0.5 

.MESH LOC=$X03 SPACING=0.2 
.MESH LOC=$X04 SPACING=0.1 
.MESH LOC=$X05 SPACING=0.05 
.MESH LOC=$X07 SPACING=0.05 
.MESH LOC=$X08 SPACING=0.05 
LOC=$X09 SPACING=0.1 
.MESH LOC=$X10 SPACING=0.25 
.MESH LOC=$X11 SPACING=0.5 
.MESH LOC=$X12 SPACING=0.5 
.-MESH LOC=$X13 SPACING=0.5 
.MESH LOC=$X14 SPACING=0.5 
.MESH LOC=$X15 SPACING=0.5 
.MESH LOC=$X16 SPACING=0.7 
.MESH LOC=$X17 SPACING=1 


mM PS PS PS PS PS DM OP OP OOOO OO OOO 
= 
ie] 
wW 
a 
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KKK KK 


-MESH LOC=0.0 
. MESH 
. MESH 
. MESH 
. MESH 


SPACING=0.05 


LOC=$Y01 SPACING=0.05 
LOC=1 SPACING=0.1 
LOC=2.25 SPACING=0.25 


LOC=$Y03 SPACING=0.5 


# Reduces mesh 


ELIMINATE 
ELIMINATE 
ELIMINATE 


density in non-critical areas to speed up simulation 


Y.DIR Y.MIN=2.25 Y.MAX=SY03 X.MIN=0 X.MAX=$X17 
Y.DIR Y.MIN=0.5 Y.MAX=2 X.MIN=0 X.MAX=14 
Y.DIR Y.MIN=2.5 Y.MAX=3.75 X.MIN=10 X.MAX=25 


# Defines the regions in ATLAS, which are equivalent to the layers in a thyris- 


tor. 
# 
# 
# 
# 


REGION NUMBER=1 
REGION NUMBER=2 
REGION NUMBER=3 
REGION NUMBER=4 
REGION NUMBER=5 


# Defi 
ELECTR 
ELECTR 
ELECTR 





# Defi 
DOPING 
DOPING 
DOPING 
DOPING 
DOPING 


Region 
Region 
Region 
Region 


PWN 





is the N1 layer 
is the Pl layer 
is the P2 layer 
is the N2 layer 
MATERIAL=SILICON 
MATERIAL=SILICON 
MATERIAL=SILICON 
MATERIAL=SILICON 
MATERIAL=SILICON 


.MIN=0 Y.MAX=SY01 X.MIN=$X02 X.MAX=$X05 
MIN=0 Y.MAX=SY02 X.MIN=$X07 X.MAX=$X16 
.-MIN=0 Y.MAX=SY01 X.MIN=$X08 X.MAX=$X11 
M O Y.MAX=SY01 X.MIN=$X12 X.MAX=$X15 


nes the location of the contacts to certain layers 
ODE NAME= 
ODE NAME= 
ODE NAME= 








anode TOP X.MIN=$X03 X.MAX=$xX04 
cathode X.MIN=$X09 X.MAX=$X10 
gate X.MIN=$X13 X.MAX=$X14 


nes the doping profile of the layers 


UNIFORM 
UNIFORM 
UNIFORM 
UNIFORM 
UNIFORM 


N.TYPE REGION=1 CONCENTRATION=$N1_ DOPING 
P.TYPE REGION=2 CONCENTRATION=$P1 DOPING 
P.TYPE REGION=3 CONCENTRATION=$P2 DOPING 
N.TYPE REGION=4 CONCENTRATION=$N2_DOPING 
P.TYPE REGION=5 CONCENTRATION=$P1 DOPING 





# Refines the grid is areas with rapoid changes in doping with distance 
REGRID DOPING LOG RATIO=2 SMOOTH=4 X.MIN=$X04 X.MAX=SX09 Y.MIN=0 Y.MAX=S$Y02 


SOLVE 


INIT 





# Saves and plots the structure file 
save outf=A 3 2.str 
tonyplot A 3 2.str 


# Defines the mobility and recombination models used for this simulation 
MODEL BIPOLAR FERMI 


# Specifies that this simulation will be controlled by current through the an- 


ode 


CONTACT name=anode current 


# Specifies the numerical technique used 


METHOD 


NEWTON 


# Opens the output file to store data generated by the following set of solu- 


tion p 


oints 


LOG OUTFILE=A 3 2.log 


# Specifies the current levels where solutions are required 


SOLVE 
solve 
solve 
solve 


INIT 


ianode=1e-19 


ianode=1. 
ianode=1. 


2e-19 
4e-19 
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# Plots the IV 


ianode=1.6e-19 


ianode=2e-19 
ianode=4e-19 
ianode=1e-18 
ianode=1e-17 
ianode=2e-17 
ianode=4e-17 
ianode=7e-17 


ianode=1le-16 istep=0.2e-16 ifinal=le-15 name=anode 


ianode=1. 
ianode=1. 
ianode=1. 
ianode=1. 
ianode=1. 


2e-15 
2e-14 
2e-13 
2e-12 
2e-11 


tonyplot A 3 _2.log 


quit 


curve 


istep=0. 
istep=0. 
istep=0. 
istep=0. 
istep=0. 


le-15 
le-14 
le-13 
le-12 
le-11 


ifinal=le-14 
ifinal=le-13 
ifinal=le-12 
ifinal=le-11 
ifinal=le-10 
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name=anode 
name=anode 
name=anode 
name=anode 
name=anode 
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